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Abstract 

The main aim of the beam energy scan (BES) program at the Relativistic Heavy-Ion Collider 
(RHIC) is to explore the quantum chromodynamics (QCD) phase diagram. The specific physics 
goal is to search for the phase boundary and the QCD critical point. We present results from 
Au+Au collisions at various energies collected in the BES program by the Solenoidal Tracker 
At RHIC (STAR) experiment. First results on transverse momentum (pr) spectra, dN/dy, 
and average transverse mass ((mr)) for identified hadrons produced at mid-rapidity for y / 3jvjv 
= 7.7 GeV are presented. Centrality dependence of dN/dy and (pr) are also discussed and 
compared to corresponding data from other energies. In addition, first results on charged hadron 
directed (wi) and elliptic flow (1)2) for ^snn = 7.7, 11.5, and 39 GeV are presented. New results 
on event-by-event fluctuations (particle ratio, net-proton and net-charge higher moments) are 
presented for ^/snn = 39 GeV. 



1. Introduction 

The Relativistic Heavy-Ion Collider at Brookhaven National Laboratory (BNL) is built 
to study the properties of a new state of matter, called Quark Gluon Plasma (QGP). 
One of the major goals of heavy-ion collision experiments is to explore the QCD phase 
diagram. The QCD phase diagram consists mainly two phases - the QGP phase, where 
the relevant degrees of freedom are quarks and gluons, and the hadronic phase. Finite 
temperature lattice QCD calculations [1] at baryon chemical potential [is — suggest 
a cross-over above a critical temperature T c ~ 170 — 190 MeV from the hadronic to the 
QGP phase. At large /is, several QCD based calculations [2] predict the quark-hadron 
phase transition to be of the first order. The point in the QCD phase plane (T vs. /i^) 
where the first order phase transition ends is the QCD critical point. The BES program 
at RHIC [3,4] aims to search the QCD phase boundary and QCD critical point. The 
QCD phase diagram can be accessed by varying temperature T and baryonic chemical 
potential /iR . Experimentally this can be achieved by varying the colliding beam energy. 
The STAR took data in the year 2010 for the beam energies y / ijvw = 7.7 GeV, 11.5 
GeV, and 39 GeV as a first phase of the BES program. 
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Table 1 

The BES energies, corresponding /ib values [9], and total events with a minimum bias (MB) trigger 
collected by STAR during the BES run in 2010. 
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Fig. 1. (color online) Left panel: Transverse momentum spectra for charged pions at mid-rapidity 
(\y\ < 0.1) in Au+Au collisions at v / 5jvjv = 7.7 GeV. The lines are the Bose-Einstein fits to the 
distributions. Right panel: (c) dN/dy normalized by (JV P art)/2 and (d) (my) — m of 7T , plotted as a 
function of collision energy. See text for details. The errors shown are the quadratic sum of statistical 
and systematic uncertainties, except for v / sjvlv" = 7.7 GeV, which has only statistical errors. 

The results presented here are based on data taken at STAR [5] for Au+Au collisions 
at y/SNN = 7.7, 11.5, and 39 GeV in the year 2010. The main detector subsystem used for 
particle identification is the Time Projection Chamber (TPC) [6]. Particle identification is 
enhanced up to higher pt with the recent inclusion of full barrel Time Of Flight (TOF) [7] 
detector. The raw yields are extracted at low-pr using ionization energy loss (dE/dx) 
from TPC, and at higher pj- using TOF information. The identified particle results are 
presented for the mid-rapidity \y\ < 0.1 region. Directed flow results are obtained using 
the Beam Beam Counter (BBC) which provides event plane determination at forward 
pseudorapidities (3.8 < \r]\ < 5.2). Event planes used for the elliptic flow results are 
provided by both TPC (\r)\ < 1) and Forward Time Projection Chamber (FTPC) [8] (2.5 
< \rj\ < 4.2). Table 1 lists the BES energies, corresponding pb values, and total events 
collected by STAR during the BES program in the year 2010. 
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Fig. 2. Left panel (color online): n~ /n + ratio at mid-rapidity (\y\ < 0.1) for central 0—5% Au+Au 
collisions at y^sjvjv = 7.7 GeV compared to previous results from AGS [10], SPS [11], and RHIC [3,12]. 
Right Panel: (a) dN/dy normalized by (-/V par t}/2 and (b) (pt) of 7r+, plotted as a function of (JV par t). 
The errors shown are quadratic sum of statistical and systematic uncertainties, except for y'Sjvjv = 7.7 
GeV, which has only statistical errors. 



2. Results 



2.1. Transverse momentum spectra at ^Jsnn = 7.7 GeV 

Figure 1 (left panel) shows the transverse momentum spectra for tt^ in Au+Au colli- 
sions at ytSNN = 7.7 GeV. The results are shown for various collision centrality classes 
as listed in the figure. The pion spectra presented here have been corrected for the weak 
decay feed-down and muon contamination. The particle production can be characterized 
by studying the dN/dy and (%) — to for the produced hadrons, where to is the mass of 
the hadron and mj = y/ to 2 + p\ is its transverse mass. These are discussed in the next 
section. 



2.2. Energy dependence of yield, (tot), and anti-particle to particle ratio 

Figure 1 (c) shows dN/dy normalized by (A f part )/2 for 7r ± in 0-5% central Au+Au 
collisions at ^snn = 7.7 GeV and are compared to previous results at AGS [10], SPS [11], 
and RHIC [3,12]. Within errors, the yields are consistent with previous results at similar 
^/sjvjv- The 7T~ /ir + ratio at ^snn = 7.7 GeV is around 1.1 (Fig. 2 left panel). Figure 1 (d) 
shows the (tot) — to for zr* in 0-5% central Au+Au collisions at ^/sWn = 7.7 GeV. 
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Fig. 3. (color online) Charged hadron v\ vs. r)/j/b ea m- The errors shown are statistical. 

The results are also compared to previous measurements at AGS [10], SPS [11], and 
RHIC [3,12]. The results from Au+Au collisions at ^/s^n — 7.7 GeV are consistent 
with corresponding measurements at SPS energies at similar ^/snn- Both dN/dy and 
{vtlt) — m are obtained using data in the measured pr ranges and extrapolating using 
a Bose-Einstein functional form for the unmeasured pr ranges. For the present mid- 
rapidity measurements, the contribution to the yields from extrapolation to the total 
yield is about 30% for ir^. 

The (tut) — to values increase with ^/snn at lower AGS energies, stay independent of 
y/ snn at the SPS and RHIC 7.7 GeV collisions, then tend to somewhat rise further with 
increasing ^/snn at the higher beam energies at RHIC. For a thermodynamic system, 
{vtlt) — to can be an approximate representation of the temperature of the system, and 
dN/dy oc ln(y / s]v r /v) may represent the entropy. In such a scenario, these observations 
could reflect the characteristic signature of a first order phase transition, as proposed 
by Van Hove [13]. However, there could be several other effects to which (tot) — to is 
sensitive, which also need to be understood for proper interpretation of the data [14]. 

2.3. Centrality dependence of dN/dy and (pt) 

Figure 2 (right panel) shows (a) the comparison of collision centrality dependence of 
dN/dy of 7r + normalized by (A par t)/2, between new results at ^/snn = 7.7 GeV and pre- 
viously published results at ^/snn = 62.4 and 200 GeV from the STAR experiment [12]. 
The yields of charged pions decrease with decreasing collision energy. At low energy, 
dN/dy/[0.5(N pliIt )} for charged pions is almost constant as a function of collision cen- 
trality. This supports the idea that particle production is dominated by soft processes at 
7.7 GeV. Right panel (b) of Fig. 2 shows the comparison of (px) as a function of (N part ) 
for 7T + from Au+Au collisions at y'SNN = 7.7 GeV with the same from collisions at 
^snn = 62.4 and 200 GeV [12]. For the collision centralities studied, the dependencies 
of (pt) on (A part ) at y/s^N — 7.7 GeV are similar to those at ^/snn = 62.4 and 200 
GeV. The values of (pt) increase from peripheral to central collisions. This indicates that 
collectivity increases with collision centrality. The (pt) values also increase with collision 
energy. 
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Fig. 4. (color online) Left panel: (a) V2 obtained from event plane method plotted as a function of pt for 
charged hadrons. The error bars include only statistical uncertainties. Right panel: (b) V2 obtained using 
the 4-particle cumulant method plotted as a function of pt for charged hadrons. The line represents the 
fit to 200 GeV data, (c) Ratio of the fit in (b) to V2 from all the energies as a function of px- 



2.4. Azimuthal anisotropy 



There are two types of azimuthal anisotropy that are widely studied in heavy-ion 
collisions, directed flow v\ and elliptic flow V2- Directed flow measurements at forward 
rapidities describe the "side-splash" motion of the collision products. The dependence 
of V\ on rj around mid-rapidity is discussed in the literature as a possible signature of a 
first order phase transition [15]. Elliptic flow provides the possibility to gain information 
about the degree of thermalization of the hot, dense medium. 

Figure 3 shows charged hadron v\ results for the 0-80% Au+Au collisions at ^snn = 7.7 
GeV, 11.5 GeV, and 39 GeV. These are compared to corresponding results from 30- 
60% Au+Au collisions at ^/sWn = 62.4, 200 GeV [16], and 0-60% Cu+Cu collisions at 
sj snn — 22.4 GeV [17]. The mid-rapidity region corresponds to the produced particles, 
while the forward rapidity corresponds to the transported particles. At mid-rapidity, 
all the results have comparable values. At forward rapidity, the trend of v\ is energy 
dependent [3]. When |ry| is scaled with the corresponding beam rapidities (j/bcam) for 
different energies, a universal curve is obtained as shown in the figure for the measured 
M/ybcam < 0.5 range. The y beam for ^/s^n = 7.7, 11.5, 39, 22.4, 62.4, and 200 GeV are 
2.1, 2.5, 3.2, 3.7, 4.2, and 5.4 respectively. 

Figure 4 (a) shows i>2 obtained from event plane method [18] plotted as a function of 
Pt for charged hadrons in 40-50% Au+Au collisions at ^/snn = 7.7, 11.5, and 39 GeV. 
The results are compared to corresponding results from y / s] v r /v = 62.4 and 200 GeV [19]. 
The figure shows that V2(pr) at 7.7 GeV is smaller than that at 11.5 GeV which is 
smaller than that at 39 GeV, suggesting that elliptic flow increases as the BES energies 
increase. Right panel (b) shows V2 obtained using the 4-particle cumulant method [20] 
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Fig. 5. (color online) Left panel: K/tt fluctuations expressed as fXdym plotted as a function of collision 
energy. The errors are quadratic sum of systematic and statistical uncertainties except for y/sjviv = 39 
GeV, which has only statistical errors. Right panel: same as (a) but for p/n fluctuations. 

plotted as a function of pt for charged hadrons in 20-30% Au+Au collisions at ^/JWn = 
39, 62.4, and 200 GeV, from STAR. The results are compared with corresponding results 
in Pb+Pb collisions at y/sWjf = 2.76 TeV, from the ALICE experiment [21]. The line 
represents the fifth-order polynomial fit to 200 GeV data. Panel (c) shows the ratio of 
the fit to V2 from all the energies as a function of px- It is observed that «2{4}(pt) for all 
the energies show similar values (within ~ 10% from 200 GeV data) beyond pr = 500 
MeV/c. This saturation is very interesting considering the wide energy range 39 GeV to 
2.76 TeV. 



2.5. Fluctuations 



The non-monotonic behavior of fluctuations in the particle ratios such as K/tt and 
p/n as a function of collision energy could indicate the presence of a QCD critical point 
or the phase transition. These ratio fluctuations are quantified by the variable fd yn [22]. 
Earlier measurements [23] of particle ratio fluctuations used the variable given by: crd yn = 

si g n (°data - cr mixod)\/l cr data- cr mixedl: where CT is the relative width of the K / ir and p/ir 
distribution in either real data or mixed events. It has been shown that Cdyn ~ ^dyn- 

Figure 5 (left panel) shows the energy dependence of dynamical K/tt fluctuations 
expressed as Odyn [24]. The 0-3.5% central Pb+Pb collisions NA49 results (solid squares) 
show decrease in dynamical fluctuations as a function of collision energy. The 0-5% central 
Au+Au collisions STAR results (solid stars) for ^siviv = 19.6, 62.4, 130 and 200 GeV, are 
constant as a function of beam energy. The newly measured K/tt fluctuation results from 
0-5% central Au+Au collisions at ^/snn — 39 GeV are similar to other STAR results [22]. 
The UrQMD model calculations (open stars) using the STAR detector acceptance, give 
approximately 4% K/tt fluctuations for all the beam energies. The right panel of Fig. 5 
shows the energy dependence of dynamical pjix fluctuations expressed as <7dyn [24]. The 
STAR results (solid stars) from 0-5% central Au+Au collisions are compared with those 
from 0-3.5% central Pb+Pb collisions from the NA49 experiment (solid squares). The new 
results from 0-5% central Au+Au collisions at ^Jsnn = 7.7 and 39 GeV are also shown. 
All the data points show increase in p/ir fluctuations as a function of beam energy. The 
STAR results from 7.7 GeV are in close agreement with those from the NA49 experiment 
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Fig. 6. (color online) Left panel: y'Sjvjv dependence of recr 2 for net-proton distributions measured at 
RHIC. Errors are quadratic sum of systematic and statistical uncertainties, except for y/sjvjv = 39 GeV, 
which has only statistical errors. Right panel: no 2 for net-charge distributions as a function of (-/V par t). 
Results from y/Sjvjv = 39 GeV are compared with those from -y/sjvjv = 200 GeV [27]. Also shown are 
the calculations from the HIJING and HRG models. 

at similar ^/smn- The UrQMD model calculations (open stars) after correcting for the 
STAR detector acceptance reproduce the increasing trend seen in data as a function of 
beam energy. The STAR results presented here are from the similar px ranges (pion, 
kaon: 0.2-0.6 GeV/c and proton: 0.4-1.0 GeV/c) for both TPC and TOF. 

It has been shown that higher moments of distributions of conserved quantities (net- 
baryon number, net-strangeness and net-charge), measuring deviations from a Gaussian, 
are sensitive to the critical point fluctuations [25]. The moments, standard deviation a, 
skewness S, and kurtosis k, of conserved quantities distributions, are defined as: a = 
y/((N-{N)f), S = ((N - (N)) 3 )/a 3 , and k = ((N - (TV)) 4 ) /a 4 - 3, respectively. 
The products of the moments such as na 2 and Sa, are related to the ratio of conserved 
quantities number susceptibilities (x) at a given temperature (T) computed in QCD 
models as : Sa ~ X^ / X^ an d Ka 2 ~ X^ / X^ ■ Close to a critical point, models predict 
the conserved quantities number distributions to be non-Gaussian and susceptibilities to 
diverge causing Sa and na 2 to deviate from constants and have large values. 

Figure 6 (left panel) shows the energy dependence of na 2 for net-proton distributions, 
compared to several model calculations that do not include a critical point. The na 2 as a 
function of the beam energy does not show any non-monotonic behavior and is consistent 
with unity. The new STAR result for na 2 from ^Snn = 39 GeV is consistent with previ- 
ous STAR measurements [26] . The right panel shows the Ka 2 for net-charge distributions 
as a function of (N part ). STAR results for ^/snn — 39 and 200 GeV are consistent with 
each other. HIJING results are consistent with the data except for the most central col- 
lisions. The na 2 for central collisions is consistent with the hadron resonance gas (HRG) 
model [28] which assumes thermal equilibrium. 



3. Summary 



The RHIC beam energy scan program has started. The aim of the BES program is to 
search for the QCD phase boundary and QCD critical point. In the first phase of the BES 
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program, the STAR experiment took data for the beam energies ^/snn = 7.7, 11.5, and 
39 GcV during 2010. The event statistics goals for three beam energies were achieved. 
The first results from the BES program are presented in this paper. Measurements of 
identified particle production at ^Snn = 7.7 GeV suggest that particle production scales 
with (-/Vp ar t). The measurements of (pt) as a function of collision centrality indicate 
that the collectivity increases with collision centrality. The mid-rapidity ir~ /tt + ratio at 
Y / s]v7V = 7.7 GeV is close to 1.1. Directed flow results for all the beam energies scale with 
??/ybcam (for measured \i]\/ybcnm < 0.5 range), extending the already established scaling 
behavior down to 7.7 GcV. The t>2{4} (pt) shows saturation above pr = 500 MeV/c 
for all the beam energies from 39 GcV through 2.76 TeV. Particle ratio fluctuations for 
the energies presented are consistent with the established trends. The net-proton and 
net-charge results are consistent with the HRG model. 
We wish to acknowledge the support from DOE. 
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